Phase-space characterization of complexity in quantum many-body dynamics 
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We propose a phase-space Wigner harmonics entropy measure for many-body quantum dynamical 
complexity. This measure, which reduces to the well known measure of complexity in classical 
systems and which is valid for both pure and mixed states in single-particle and many-body systems, 
takes into account the combined role of chaos and entanglement in the realm of quantum mechanics. 
The effectiveness of the measure is illustrated in the example of the Ising chain in a homogeneous 
tilted magnetic field. We provide numerical evidence that the multipartite entanglement generation 
leads to a linear increase of entropy until saturation in both integrable and chaotic regimes, so 
that in both cases the number of harmonics of the Wigner function grows exponentially with time. 
The entropy growth rate can be used to detect quantum phase transitions. The proposed entropy 
measure can also distinguish between integrable and chaotic many-body dynamics by means of the 
size of long term fluctuations which become smaller when quantum chaos sets in. 
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I. INTRODUCTION 

Understanding the dynamics of quantum systems is 
a challenging task of immense importance in a variety 
of fields including condensed matter physics and quan- 
tum information science. Quantum dynamical complex- 
ity refers to the lack of a simple description of the evo- 
lution of a quantum system. From a computational per- 
spective, it implies the inevitable loss of predictability 
of system evolution using classical simulation. In many- 
body interacting quantum systems, complexity can be 
attributed to non-integrability or to the tensor product 
structure of the Hilbert space. Hence, quantum chaos 
and entanglement have deep implications in characteriz- 
ing quantum many-body dynamical complexity. 

In classical physics, it is very well known that there 
exists a direct correlation between chaos and complex- 
ity. Classically chaotic systems are characterized by ex- 
ponentially diverging nearby trajectories, with a rate de- 
termined by the Lyapunov exponent. Complexity then 
arises from the fact that the orbits of such deterministic 
systems are completely random and unpredictable with 
positive algorithmic complexity [l| . In quantum mechan- 
ics, trajectories in standard treatments are forbidden by 
the Heisenberg uncertainty principle and therefore the 
above notion of complexity cannot be directly translated 
to quantum systems. 

However, the phase-space approach can be equally 
used for both classical and quantum mechanics. In the 
context of classical systems, it has been shown that the 
second moment of the Fourier components of the classi- 
cal distribution function grows linearly for an integrable 
system while it grows exponentially for a chaotic sys- 



tem, with a rate determined by the Lyapunov exponent 
characterizing the local exponential instability. Thus, the 
growth rate of the second moment of Fourier components 
(harmonics) is a good measure of the complexity of clas- 
sical dynamics [2j. In a similar way, for single-particle 
quantum systems the second moment of harmonics of the 
Wigner distribution function of a quantum state^pure or 
mixed, is a measure of quantum complexity [2|, [3|. Note 
that in quantum systems with few degrees of freedom an 
exponential growth of the number of harmonics is possi- 
ble only up to the Ehrenfest time scale, after which the 
growth is at most linear [3]. Moreover, the number of 
harmonics of the Wigner function can be used to detect, 
in the time domain, the crossover from integrability to 
chaos [1]. 

For quantum many-body systems the situation is more 
complicated. First note that quantum dynamical en- 
tropies, which generalize the Komologrov-Sinai entropy 
to quantum dynamical systems, can be positive even for 
integrable dynamics This behavior may appear, at 
least at first sight, somehow surprising since in classi- 
cal dynamics positive Komologorov-Sinai entropy implies 
chaos. Another interesting property is that, as shown 
in Ref. [6;], the rank of the matrix product operator 
representation of the pure quantum states in the time- 
dependent density-matrix rcnormalization group, typi- 
cally grows exponentially even for integrable system with 
finite number of particles. This inefficiency of the clas- 
sical simulation of many-body quantum dynamics can 
be attributed to entanglement and is consistent with the 
linear growth of the entanglement block entropy for in- 
tegrable spin chains Q. 

Several very interesting definitions of quantum com- 
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plexity have been proposed, e.g., see Ref. @ and ref- 
erences therein. On the other hand, to the best of our 
knowledge none of them satisfies all the following require- 
ments, which a notion of complexity should possess in 
order to be both meaningful and practically useful: 

(i) To provide a unified description of both one- and 
many-body dynamics; 

(ii) To reproduce at the classical limit the well-known 
notion of classical complexity based on the local 
exponential instability of chaotic dynamics; 

(iii) To be applicable to both pure and mixed states; 

(iv) To be practically useful, that is, convenient for nu- 
merical investigations. 

The purpose of the present paper is to propose a notion 
of complexity that fulfills the above criteria. By extend- 
ing previous investigations [2HJ| to many-body quantum 
dynamics, we propose the number of harmonics of the 
Wigner function as a suitable measure of complexity of 
a quantum state. Indeed, as the phase-space formula- 
tion of quantum dynamics can be directly generalized to 
many-body systems, the harmonics of the Wigner func- 
tion seem to be very promising in quantifying the com- 
plexity of many-body quantum systems as well. Hence, 
in this paper, we introduce a Wigner harmonics entropy 
measure of complexity and then illustrate its usefulness 
by means of numerical simulations carried on a paradig- 
matic spin-chain model, the Ising chain in a tilted mag- 
netic field. We will show that the entropy grows linearly 
until saturation in both integrable and chaotic regimes, 
so that in both cases the number of harmonics of the 
Wigner function grows exponentially with time. We will 
provide numerical evidence that this growth must be at- 
tributed to multipartite entanglement generation. Our 
results demonstrate that the growth rate can be used 
also to detect quantum phase transitions. Finally, the 
proposed entropy measure can also distinguish between 
integrable and chaotic many-body dynamics, by means of 
the size of long term fluctuations, which become smaller 
when a transition to chaos occurs. 

This paper is organized as follows. In Sec. [Ill we define 
our phase-space quantum complexity measure, based on 
the harmonics of the Wigner function. The working of 
such measure is illustrated in the dynamics of a many- 
body spin-chain model, introduced in Sec. IHII and inves- 
tigated in detail in Sec. IIV1 Finally, our conclusions are 
drawn in Sec. IV] 



II. HARMONICS OF THE WIGNER FUNCTION 

The phase-space representation of quantum mechan- 
ics is a very enlightening approach as it allows a direct 
comparison between quantum and classical dynamics. In 
particular, the complexity of a quantum state or of a 



classical distribution function can be measured by the 
richness of their phase space structure. 

In the quantum case, the phase-space approach to 
complexity is particularly convenient for systems whose 
Hamiltonian can be written in terms of a set of bosonic 
creation-annihilation operators: 



H(a\, ... ,0^,0,1, ...,a N ;t) = IT (0) (ni, h N ) 
+ H ( - 1 \a\, ajy, ax,..., a N ;t), 



(1) 



with [&i, a,j] — [a|, at] = 0, [a|, aj] = dy, and the number 

operators hi = fit a,. 

We will use the method of c-number a-phase space 
borrowed from quantum optics (see for example Ref. [8j]). 
The Wigner function W(a, a*;t) of a state p(t) is defined 
by 



W{a,a*;t) 



1 



ir 2N h N 



d rjexp 



(2) 



where rj — (rji, t]n) and a. = (a±, ...,apf) are TV- 
dimensional complex variables, the integration runs over 
the complex 7/i-planes for i — l,...,iV, the displacement 
operator 



D (rj) = exp 
and the coherent states 



' N 
.i=l 



\a) = |aia2-.Q!jv) = D 



I00....0), 



(3) 



(4) 



with |a!j) being eigenstate of the annihilation operator 
Oj, i.e., ai\a>i) = -^=\ai), and |00...0) being the vacuum 

state. We define the harmonic's amplitudes W m (l;t) of 
the Wigner function by the iV-dimensional Fourier ex- 
pansion 



(5) 



where m, I, 6 are A^-dimensional vectors, whose compo- 
nents Ik > 0, < 9k < 27r are defined by the relations 
otk = VTk~e ~ l8k , k = 1,...,N. Here Ik and 9k can be re- 
garded as our quantum phase space variables, analogous 
to the action and angle variables in the classical phase 
space. Note that W- m = W^. The Wigner function's 
normalization condition J d 2 otW(a, a* ; t) — 1 simply 
implies that J dIWo(l;t) = 1, while there are no restric- 
tions on when m/0. 

In Refs. |3|, |4[ , the number of harmonics of the Wigner 
function was estimated by yj (m 2 ) t , with (m 2 ) t being the 
second moment of the harmonics distribution: 



(m 2 ) t 



^m 2 W m (t), 



(6) 
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where 



III. THE MODEL 



W m (i) 
The harmonics 

E m Wm = 1. 



JdI\W m (I;t)f 
E m Jdl\W m (I;tW 



(7) 



distribution W m is normalized, 
For one-body systems, the second 



moment (m ) t provides a reliable estimate of the num- 
ber of harmonics in a generic chaotic case [3j and is able 
to distinguish, in the semiclassical region, between intc- 
grable and chaotic regimes Q. In the first case, ■J (m 2 ) t 
grows linearly in time, in the latter exponentially. On 
the other hand, we expect that the number of harmonics 
always captures the complexity of motion, including the 
case of many-body systems without classical analogue. 
Hence, we propose as a complexity measure the entropy 



S(t) 



W m (t)HW m (t)] 



(8) 



where the sum over m is limited to mi, mjv > since 
harmonics W m and W- m are not independent but triv- 
ially related by the relation W m = W_ m (note that the 
same limitation must now be taken in Eq. (J7J in order 
to properly normalize the distribution W m ). The num- 
ber of harmonics of a generic state p{t) can therefore be 
measured by exp[S(£)]. For the models discussed in this 
paper both the second moment (m 2 ) and the entropy S 
provide qualitatively the same results @. 

The main computational advantage of the above c- 
number a-phase space approach is that the Wigner func- 
tion's harmonics W m can be computed very conveniently 
from the density matrix written in the basis of the eigen- 
vectors |n) = |ni...7ijv) of the unperturbed Hamiltonian 
H(°\ Indeed, using the well-known matrix elements of 
the displacement operator fioj . 



(rii+mi\D(r)i)\ni 



(9) 

{m,mi > 0, i = 1, ...,JV), where L™*(x) is a Laguerre 
polynomial, the ^-integration in Eq. ^ can be car- 
ried out explicitly. After that, using the orthogonal- 
ity and completeness properties of the Laguerre poly- 
nomials along the lines of Ref. Q, we can express the 
Wigner harmonics W m (l;t) in terms of the matrix ele- 
ments (n + m|/5(i)|n) and finally obtain 



W m (t) 



E n |<n + m|p(t)|n) 



£ n |(n + m|p(t)|n) 



(10) 



Finally, we point out that our approach remains valid 
also for classical systems, provided the Wigner function 
is substituted by the classical phase-space distribution 
function in the a^-coordinates, with at = \/Tk~e~' iek , 
{Ik, 9k} being a set of action-angle variables for the un- 
perturbed, integrable Hamiltonian H . 



In order to investigate the working of our complex- 
ity measure, we consider, as an illustrative example, the 
Ising chain of N spins in a tilted magnetic field. The 
Hamiltonian reads 



H 



(11) 



where J is the spin-spin coupling constant, af are the 
Pauli operators for the i-th. spin, and h x , h z are the field 
amplitudes along x and z directions, respectively. We 
set ft, = J = 1. This chain is in general non- integrable, 
except for the two integrable limits h x = or h z = 0. 
The integrable model h z = corresponds to the Ising 
model in a transverse field and exhibits a quantum phase 
transition at J = h x [lll |. 

Using the Schwinger boson representation [121 ], the 
above spin Hamiltonian is mapped onto an interacting 
boson Hamiltonian. Each spin operator bi at the site i is 
replaced by two Schwinger bosons, fij and bi, correspond- 
ing to spin up {t} and down {!}. The spin operators can 
be represented as follows: 



it; 



erf = d]6j, and a,- = Ma,, 



t; 



(12) 



where a 



\{b x ± ia y ). Since we have spin- 1/2 par- 
ticles, the physical subspace is singled out by the con- 
straints n a i + nt,i = 1 (i = l,...,iV), where n a i and 
nbi denote the number of up and down spins at site i 
(n a i = a\a,i, hbi = b\bi)- Now, Eq. (fTTj) takes the form 



N-l 



i=l 
N 

+ Y}h x {a% + bla,) + h z {a\a t - &]£;)]. (13) 

i=l 

As the Hamiltonian is now expressed in terms of a set of 
bosonic creation-annihilation operators, it then follows 
that the above explained phase-space approach can be 
used to probe the dynamical complexity of the spin chain. 

For a chain of N spins, n and m in Eq. (ITU1) are 2N- 
dimensional vectors, 

n = (n a i,n b i,n a2 ,n b2 , ...n aN ,n bN ), 
m = (m a i,m bl ,m a2 ,m b2 ,...m aN ,m bN ), (14) 

where the first subscript refers to the spin type and sec- 
ond refers to the spin site. The possible values of n's are 
and 1 with the constraint n aj ; + n b i = 1. Similarly, the 
possible values m's can take are —1,0 and 1 with the con- 
straint m a i + m b i = to remain on the physical subspace. 
Indeed, the following cases are possible: (i) m a i = — 1, 
m b i — 1, corresponding to the transition of the i-th spin 
from up to down, (ii) m a i = 1, m b i — — 1 (transition of 
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the i-th spin from down to up), and (iii) m ai — mu — 
(no transition for the i-th spin). Due to the trivial rela- 
tion W-m = W m , we limit the summation ([5]) to m ai > 
(and, consequently, < 0), for i = 1, N, in order to 
take into account independent terms only. That is, for 
a chain of N spins, only 2 N values of the W m are inde- 
pendent and are considered for calculations. Note that 
the maximum possible value of the entropy measure S(t) 
is iVln2. This value is reached when maximum mixing 
occurs so that all the harmonics are equally distributed. 



IV. PHASE-SPACE CHARACTERIZATION OF 
COMPLEXITY 

A. Initial Growth of S(t) 

In this section, we study in detail the time evolution 
of the entropy S(t) at small times. The initial state is 
chosen to be a pure state with all spins pointing down- 
ward in the z direction, i.e., |^i n ) = | W ■■■ -l)- From the 
definition of Wigner harmonics W m (i) in Eq. (p~0|) . it is 
clear that only the m= (0,0, ...0) harmonics component 
is excited, the initial number of harmonics is equal to 
unity and hence S(t = 0) = 0. 

To relate the growth rate of the harmonics to the com- 
plexity of the dynamics as in Ref. Q, first we consider 
the short time behavior of the entropy measure S(t) when 
the system undergoes a transition to quantum chaos as 
detected by a change of Poisson to Wigner distribution 
in the level statistics. To that end, let us rewrite the 
Hamiltonian H as in Eq. ([1]) with 

i i i 

Here is the integrable Hamiltonian and H W repre- 
sents the perturbation to the chain induced by an exter- 
nal transverse field. As the perturbation is increased, a 
transition to Wigner-type level statistics and hence quan- 
tum chaos occurs. In particular, for h x = h z = J the 
system can be considered as fully chaotic [6J. 

The initial state is an eigenstate of the unperturbed 
integrable Hamiltonian H^°\ Hence with zero perturba- 
tion (h x = 0), there is no evolution and S(t) remains zero 
at all times. However, by adding a small perturbation 
to the system i.e., a small transverse field h x , the initial 
state is no longer an eigenstate of the Hamiltonian H and 
transitions to many other states occur. Then, besides the 
zeroth harmonics i.e., m = (0,0, ...0), higher harmonics 
are also excited and the entropy S(t) increases. Smaller 
the transverse field, less complex is the dynamical evo- 
lution and therefore a lower value of the growth rate for 
entropy S(t) is obtained. This is clearly seen from our 
results in Fig. [Ha). For instance, with h x = 0.2, the 
chain is near the integrable regime and S(t) is about 1.0 
at time t = 1.0. By contrast, in the chaotic regime with 
h x = 0.9, S(t) increases to 6.3. Interestingly, for cases 




FIG. 1: (color online) Time dependence of the entropy mea- 
sure S(t) for a chain of 10 spins with (a) non-integrable Hamil- 
tonian as given in Eq. (|15[) with a longitudinal field h z = 1.0 
and (b) integrable Hamiltonian as given in Eq. (|16|) . Curves 
from top to bottom correspond to transverse field h x = 1.0 
to 0.1 in decreasing steps of 0.1. Note that during an ini- 
tial time window S(t) is clearly seen to grow linearly with 
time, implying the exponential growth of the number of har- 
monics in both (a) and (b), for a sufficiently large h x . All 
the parameters mentioned here, and in the other figures, are 
dimensionless (we set % = J = 1). 

with a sufficiently strong perturbation, S(t) is seen to in- 
crease linearly with time within a time window (t < 0.6). 
Such linear increase of S(t) implies an exponential growth 
of the number of harmonics. 

Consider then another situation for the Ising chain, 
with the magnetic field applied in the transverse direction 
only: 

# (0 WX;5f5? +1 ; H^=J2 h ^f- ( 16 ) 

i i 

Note that this Hamiltonian is integrable for all values of 
the perturbation h x [l3j . On the basis of previous find- 
ings in few-body problems @, 0| ; we might expect a linear 
increase of the number of harmonics, corresponding to a 
logarithmic growth of S(t). On the contrary, as shown 
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integrable and integrable spin chains. We must therefore 
seek an underlying mechanism to account for this some- 
what counter-intuitive behavior of many-body quantum 
systems. 



B. Wigner Harmonics and Entanglement 

A source of quantum complexity in many-body sys- 
tems is the entanglement due to the interaction between 
the different constituent parts. Recent studies indicated 
that some measure of the entanglement entropy can also 
grow linearly with time Q. We will therefore inquire 
whether or not the lack of distinction between integrable 
and non-integrable models shown above is related to the 
generation of multipartite entanglement. 

To quantify the extent of multipartite entanglement 
generated in a spin chain, we adopt the multipartite en- 
tanglement measure used in Ref. [l4|. Specifically, the 
system under consideration is partitioned into two sub- 
systems A and B, made up of ua and Ub spins, respec- 
tively. The participation number Nab, defined as the 
reciprocal of the purity of one of the two subsystems, 
i.e., 



N A i 



(17) 



FIG. 2: (color online) Dependence of entropy S(t) at time 
t = 1 on the strength of the external perturbation h x for 
different chain length N for (a) a non-integrable spin chain 
and (b) an integrable spin chain. For both cases and for 
sufficiently large h x , it is seen that S(t = 1) scales linearly 
with N. 



in Fig. QJb) the short time behavior of S(t) is linear. 
A comparison between Fig. Ufa) and Fig. QJb) clearly 
shows that one cannot distinguish between a chaotic and 
an integrable many-body system by solely examining the 
initial growth of the number of Wigner harmonics. 

To gain more insights, we next study how the initial 
linear growth rate of S(t) depends on the number N of 
spins in the chain, for both the non-integrable model (|15[) 
and the integrable model (I16|) . In particular, we vary N 
from N — 8 to N — 14. Figure [U[a) shows the value 
of S(t) at a fixed time t = 1 as a function of h x , for 
the non-integrable model with four different values of 
N. Note that S(t = 1) can be understood as the av- 
erage entropy production rate for t £ [0,1]. It can be 
observed that S(t = 1) scales with N linearly. For ex- 
ample, for h x = 1.0, S(t — 1) increases by a constant 
value (s=s 1) as N increases in steps of two. Interestingly, 
as shown in Fig. HJb), exactly the same behavior is ob- 
served for the integrable model. This further strengthens 
our early finding that the initial growth of the number 
of Wigner harmonics is qualitatively the same for non- 



accounts for the bipartite entanglement between A and 
B. Here, pa is the reduced density matrix of subsystem 
A. The physical meaning of Nab is that it effectively 
counts the relevant terms in the Schmidt decomposition 
of the total wavefunction into the sum of direct products 
of wavefunctions of the two subsystems. The mean value 
of Nab, averaged over all possible partitions, quantifies 
the degree of multipartite entanglement in the system, 
while its variance measures how well the entanglement is 
distributed. For systems of large size N » 1, the sta- 
tistical weight of unbalanced partitions becomes negligi- 
ble [141 ] and hence only balanced partitions are considered 
here. 

In Fig. [3] we present the time dependence of the mean 
value of the participation number (Nab), starting from 
the same initial state as before, i.e., |\^m) = I 11 ■■■ I), 
for both integrable and non-integrable models. The ini- 
tial state is not entangled and hence (NAB)(t = 0) is 
given by its minimum value: unity. Entanglement is then 
generated by the dynamical evolution of the spin chain 
and hence (Nab) increases. Remarkably, after a short 
time interval (t < 0.4) and for a sufficiently large value 
of h x , {Nab) reaches the saturation value almost expo- 
nentially fast, for both the integrable and non-integrable 
models. Though the production of entanglement is some- 
what slower for the integrable chain as compared to the 
non-integrable case, an exponential-like fast increase of 
{Nab) is seen in both situations. To visualize this more 
clearly, we plot an exponential fit for h x — 0.8 in both 
panels of Fig. [3] These results indicate that we can as- 
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FIG. 3: (color online) Time dependence of average value of 
participation number {Nab} calculated over all balanced bi- 
partitions of the system for (a) non-integrable and (b) inte- 
grable model with parameters discussed in Fig. [1] Curves 
from top to bottom correspond to transverse field h x = 1.0 
to 0.1 in decreasing steps of 0.1. Within a small time win- 
dow, (Nab) oc e At for relatively large h x , as shown by an 
exponential fit (circles) for h x — 0.8 in both panels. 



FIG. 4: (color online) Comparison of the dynamics of normal- 
ized entropy S'norm and global entanglement GE for a chain 
of 10 spins with transverse field h x = 0.8. Panel (a) cor- 
responds to the non-integrable model with longitudinal field 
h z = 1.0 and panel (b) corresponds to the integrable model. 
Here S'norm is the entropy S divided by its maximum value 
iVln2. A close correspondence between the dynamics of two 
measures is evident. 



cribe the exponential growth of Wigner harmonics to the 
fast entanglement generation in the chain. 

To better clarify this latter point, we have compared 
the time-dependence of S(t) with the so-called "global 
entanglement" (denoted GE) [l5[. In Ref. [l6j], it was 
shown that GE is related to the averaged one-qubit pu- 
rity, i.e., 

GS = 2^1-l^Tr[^]V (18) 

where pk is the density matrix of the fc-th spin after trac- 
ing over all other spins in the system. GE is the average 
bipartite entanglement over all possible bipartitions be- 
tween a single qubit and the rest of the system. It is 
easy to see that < GE < 1. Values of GE close to 
one indicate highly entangled many-body states. When 
a many-body state is not entangled, GE equals to zero. 

For the initial state l^m) = | 11 ■■■ I), we present in 
Fig. |4]a comparison between S and GE. To better visu- 



alize their similarities, we plot a normalized (to unity) en- 
tropy S'norm = S/N In 2. It is clearly observed that these 
two quantities show a high degree of resemblance in their 
time dependence. Their oscillating patterns are quite 
close and in some regimes they are almost on top of each 
other. The similarity between S'norm and GE constitutes 
strong evidence that our entropy measure S(t), though 
originated from considerations of phase-space complex- 
ity, also reflects the degree of multipartite entanglement 
in many-body systems. 

Though, in general, GE may not distinguish between 
different classes of multipartite entangled states, it is an 
indicator of the critical point of, for instance, the quan- 
tum phase transition for the Ising chain in a transverse 
magnetic field [l7[. Therefore, it is also interesting to 
investigate the behavior of S(t) in the neighborhood of 
a quantum critical point which, for the transverse Ising 
chain in Eq. (|16p with coupling strength J = 1, is at 
h x — 1. In Fig. we show the behavior of S(t = 0.5) 
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FIG. 5: (color online) Panel (a) shows S(t = 0.5) normalized 
by a factor of N In 2 as a function of the transverse field h x 
for a transverse Ising chain of 10 spins. Panel (b) shows the 
parallel results for GE at t = 2.0. Here A = h x /(J + h x ). 
Similar to GE, the complexity measure S is seen to peak 
clearly at the critical point A = 1/2. 



as well as GE(t — 2), as a function of A = h x / (J + h x ). 
Note that S and GE are plotted at different times be- 
cause of their different saturation times. In addition, in 
our calculations of S on the large-field side (A > 1/2) the 
x axis is used as the quantization axis of the basis states: 
The magnetic field term is dominant and correspondingly 
we define = Ei^f, = ^EiW+n and 

Schwinger bosons such that of = a|dj — &|5j. It is quite 
natural to consider as preferential basis the one associ- 
ated with the dominant term in the Hamitonian: the z- 
basis when A — > and the x-basis when A — > 1, and quan- 
tum phase transition corresponds to the switching from 
one preferential basis to the other. Consistent with the 
expectation that the quantum phase transition occurs at 
A = 1/2, Fig. [5ja) shows that S noim (t = 0.5), a measure 
of the growth rate of the number of Wigner harmonics, 
exhibits a sharp peak at A = 1/2. The A-dependence of 
GE shown in Fig.(5Jb) is analogous to what we observe in 
Fig-EJa). This further demonstrates the close connection 
between our complexity measure S(t) and the global en- 



tanglement GE and, in particular, the role of many-body 
entanglement in the initial growth of S(t). 

Two additional aspects of S(t) are in order. First, if 
we stick to the z axis as the quantization axis of the 
basis states, then it is found that right after the critical 
point, S(t) (if averaged over a time window to remove 
fluctuations) will show clear saturation behavior, which 
is in contrast to the monotonous increase of S(t) before 
the critical point. Second, if we switch the quantization 
axis from x to z at other values of A, then the value of 
S(t = 0.5) jumps discontinuously due to the change of 
the basis states. These additional results further suggest 
that the critical point for quantum phase transitions can 
be detected by S(t). 

Note that a different phase-space measure [IH has been 
used in the literature to detect quantum phase transi- 
tions However, such measure accounts for the extent 
at which the phase space is covered by the Husimi distri- 
bution and therefore it does not appear clear how to ex- 
tend it to a suitable complexity measure for mixed states. 
In contrast, our measure which is based on the richness 
of the phase space structure rather than on phase-space 
coverage can be used for both pure and mixed quantum 
states. For instance, it could be used also to investigate 
thermal phase transitions. 



C. Wigner Harmonics, Chaos, and Thermalization 

Our results so far indicate that due to the dynam- 
ically generated many-body entanglement, the initial 
time-dependence of S(t) does not reflect the peculiarity 
of quantum chaos in many-body quantum systems: it be- 
haves similarly in integrable and non-integrable models. 
Note that this does not contradict with previous find- 
ings regarding rapid bipartite entanglement generation in 
classically chaotic systems with two degrees of freedom 
[20| . Indeed, in systems with two degrees of freedom, the 
Hilbert space is only a product of two subspaces and the 
quantum dynamics can only generate bi-partite entan- 
glement. The rate of entanglement growth within such 
a fixed product of two subspaces is connected with the 
underlying classical dynamics. By contrast, in a many- 
body system such as our model used here the dynamics 
emanating from a local initial state is seen to explore 
more and more the tensor-product structure of the total 
Hilbert space and hence entangle more and more degrees 
of freedom during the time evolution. Our entropy mea- 
sure S(t) then is indicative of the growth rate of the num- 
ber of degrees of freedom that have been entangled by the 
dynamics (a property absent in few-body systems). 

Since the short-time behavior of S(t) is seen to be un- 
related to quantum chaos, we now examine the mani- 
festation of quantum chaos in the long-time behavior of 
S(t). This can be justified because after all, the pecu- 
liar spectral statistics of a quantum chaotic many-body 
system reflects the long-time properties of the system. 

The time-dependence of S(t) for a non-integrable 
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FIG. 6: (color online) Time evolution of entropy S(t) for a 
chain of 12 spins with (a) non-integrable Hamiltonian in Eq. 
(|15|l for h z = 1.0 and (b) integrable Hamiltonian in Eq. (|16|) . 
Curves from top to bottom corresponds to transverse field h x 
= 1.0, 0.5 and 0.2. Note that the long term dynamics of the 
entropy in the two panels is qualitatively different. 

model with the Hamiltonian given in Eq. ()15[) is shown 
in Fig. Eta). The parallel result for an integrable model 
defined by Eq. (fT6|) is shown in Fig. (6][b). In both cases, 
the time scale under study is now 20 times longer than 
that used in Fig. 1. It is seen that, in both panels (a) 
and (b), S{t) initially quickly increases and then displays 
saturation with rich oscillating behavior. The saturation 
plateau of S(t) increases as the value of the transverse 
field increases. Qualitatively, the saturation plateau can 
be attributed to an effective dimension of the Hilbert 
space that can be explored for a particular strength of 
the transverse field. To quantitatively describe this ob- 
servation, we calculate the average value of the entropy 
from time t\ = 5 to t2 = 100 as S = \ J t ^ dtS(t), where 
t = <2 — t\. The values of t\ and t% are chosen so that 
the saturation plateau is reached before t\ and ti is large 
enough to allow averaging over many oscillations of S(t) 
in the time interval between t± and ti. We have checked 
that other such choices of t\ and t 2 do not affect any 
of our observations reported below. In Fig. we plot 
S as a function of the transverse field, for the integrable 



and non-integrable models considered in Fig. [B] It is seen 
that, for small values of h x , there is a difference between 
integrable and non-integrable dynamics. However, as the 
strength of the transverse field increases, this difference 
reduces. This is somewhat expected due to the above- 
discussed many-body entanglement generation. 
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FIG. 7: (color online) Time averaged entropy, denoted S, vs 
the strength of the transverse field h x , for a chain of 12 spins. 
The non-integrable model corresponds to the Hamiltonian in 
Eq. (|15[) with h z — 1.0; the integrable model corresponds 
to the Hamiltonian in Eq. I|16p. The difference in S be- 
tween integrable and non-integrable dynamics decreases with 
increasing hx- 

In order to distinguish between integrable and non- 
integrable cases, we are thus forced to look into the os- 
cillating behavior (rather than the average behavior) of 
S(t). Indeed, from Fig. [6] one observes that in the non- 
integrable case, the oscillation amplitude of S{t) clearly 
decreases with the value of h x . The oscillation pattern 
also becomes erratic as the system gets closer to the 
chaotic regime. By contrast, in the integrable case the 
opposite trend is observed. Regular and strong quantum 
revivals in S(t) become more apparent as h x increases. 
To quantitatively describe this clear difference, we cal- 
culate the standard deviation of S(t) around the mean 
value S, 



*[S] = ^j\t[S(t)-S]2. (19) 

The results are shown in Fig. [8j With increasing pertur- 
bation, the standard deviation a[S] increases and then 
saturates in the integrable model, so that the relative size 
a[S]/S of fluctuations remains nearly constant. However, 
in the non-integrable model, the standard deviation <r[S] 
and, more markedly, a[S]/S decrease during the regular- 
to-chaotic crossover (The same qualitative behavior is ob- 
tained when the number of the spins N in the chain is 
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varied). For N = 12, a[S] in the integrable model with 
h x = 1 is around 0.878. By sharp contrast, in the non- 
integrable case with h x — 1, a[S] — 0.1495, which is 
smaller than the first case by more than 5 times. This is 
a dramatic difference considering that the total number 
of spins in the chain is only 12. 



CO 






1.0 



tion value X(t) = (a x )(t). Our results show that in the 
integrable model, the standard deviation increases with 
the transverse field whereas for the non-integrable case, 
the standard deviation decreases. This behavior, illus- 
trated in Fig. [SJ is qualitatively similar to the behavior 
of S(t) shown in Fig. [5] This shows that our complexity 
measure is related to the thermalization properties of the 
system. 
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FIG. 8: (color online) Standard deviation a[S] in the entropy 
S(t) vs the strength of the transverse field h x for a chain of 
12 spins. Here the integrable and non-integrable models are 
the same as in Fig. [7] It is observed that as h x increases, the 
standard deviation generally decreases in the non-integrable 
model, but increases in the integrable model. For large values 
of h x , <t[S] for a non-integrable chain is much smaller than 
that for an integrable chain. 

The large value of the standard deviation for the inte- 
grable model can be accounted for by the lack of ther- 
malization. Indeed, our expectation is that the onset of 
chaos leads to internal dynamical thermalization [2l[ , so 
that a statistical description is possible even though we 
have a closed, finite Hamiltonian system. Since the den- 
sity of many-body energy levels grows exponentially with 
the number of particles, even a weak interaction between 
particles typically leads to a strong mixing on nonintcr- 
acting many-body states, thus resulting in chaotic eigen- 
states. That is to say, the components of such eigenstates 
can be treated as random variables and therefore statis- 
tical methods can be applied to the description of local 
observables, in spite of the fact that close systems are 
under consideration. In such situation, fluctuations of 
the expectation values of local observables are small. On 
the contrary, in the integrable regime the lack of ther- 
malization allows large fluctuations. To verify the above 
expectations, we have considered the Pauli operator a x 
(note that, due to translational invariance of model and 
initial condition, (af ){t) is independent of i at any time 
t). Observing the long term dynamics of the chain, we 
compute in the same way as in Eq. (|19l) for S(t) the 
standard deviation a[X] of the x-polarization expecta- 



FIG. 9: (color online) Standard deviation of the x- 
polarization expectation value as a function of transverse field 
h x for a chain of 12 spins. Here the integrable and non- 
integrable models are the same as in Fig. [7] Similar to the 
case of entropy S(t), with increasing perturbation, the stan- 
dard deviation <j[X] decreases in the non-integrable model 
and increases in the integrable model. 



V. CONCLUSIONS 

In this paper, we propose an entropy measure S(t) 
for many-body quantum-dynamical complexity, by ex- 
tending the Wigner harmonics measure introduced in 0- 
4] for single-particle quantum dynamics. The effective- 
ness of this measure is illustrated in the example of the 
Ising chain in a homogeneous tilted magnetic field. The 
Wigner harmonics entropy S(t) exhibits an initial lin- 
ear growth in both integrable and chaotic regimes, until 
saturation occurs due to the finite size of the Hilbert 
space. Therefore, in both integrable and chaotic regimes 
the number of harmonics of the Wigner function grows 
exponentially with time. In classical dynamics, an expo- 
nential growth of the number of harmonics of the classi- 
cal phase-space distribution function implies chaotic dy- 
namics. Therefore, the observed exponential growth of 
Wigner harmonics in the many-body quantum integrable 
regime must be attributed to a source of complexity ab- 
sent in classical dynamics, that is, entanglement. We 
have numerically demonstrated the close connection be- 
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twccn our complexity measure S(t) and multipartite en- 
tanglement, thus providing evidence that the initial lin- 
ear growth of S(t) has to be ascribed to multipartite 
entanglement generation. The Wigner harmonics mea- 
sure S(t) can also distinguish between integrable and 
chaotic many-body systems, by means of the size of long 
term fluctuations, which are smaller in the chaotic regime 
where a statistical description of the system is legitimate 
and the relative size of fluctuations drops when the sys- 
tem size increases. 

The main advantage of the phase-space approach to 
complexity resides in its generality. At the classical limit, 
the harmonics of the phase-space distribution function 
reproduce the well-known notion of complexity based on 
local exponential instability [|| : the number of harmonics 
grows linearly for integrable systems and exponentially 
for chaotic systems. In single-particle quantum mechan- 
ics, an exponential growth of the number of harmonics is 
possible only up to the Ehrenfest time scale, after which 
the growth is at most linear Q. Furthermore, the num- 



ber of harmonics of the Wigner function can be used to 
detect, in the time domain, the crossover from integra- 
bility to chaos 4]. For quantum many-body systems, 
the Wigner harmonics entropy measure S(t) proposed 
in this paper signals the generation of multipartite en- 
tanglement and can be used to detect quantum phase 
transitions. In relation to other measures of complexity 
based on the efficiency of the best classical simulations 
of quantum systems [£jj , our approach has the advantage 
that it does not rely on a specific computational method 
like the time-dependent density-matrix renormalization 
group. Finally, we point out that, in contrast to other 
quantum phase-space a ppr oaches based on the moments 
of the Husimi function [l8l . fl9| , our complexity measure 
works equally well for either pure or mixed quantum 
states. Therefore, our measure could be studied in re- 
lation to mixed-state entanglement. This would be par- 
ticularly interesting as mixed-state entanglement is at 
present not well understood and is the focus of ongoing 
research. 
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